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Abstract 

We explore the Higgs and sparticle spectroscopy of supersymmetric S'f/(4)c x 
SU{2)l X SU{2)ji models in which the three MSSM gauge couplings and third fam- 
ily (t-b-r) Yukawa couplings are all unified at Mqut- This class of models can be 
obtained via compactification of a higher dimensional theory. Allowing for opposite 
sign gaugino masses and varying rrit within la of its current central value yields a 
variety of gauge- Yukawa unification as well as WMAP compatible neutralino dark 
matter solutions. They include mixed bino-Higgsino dark matter, stau and gluino 
coannihilation scenarios, and the A-resonance solution. 
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1 Introduction 

Supersymmetric (SUSY) 5*0(10) GUT (grand unified theory), in contrast to its non- 
SUSY version, yields third family (t-b-r) Yukawa unification via the unique renormal- 
izable Yukawa coupling 16 ■ 16 ■ 10, where the 10-plet is assumed to contain the two 
minimal supersymmetric standard model (MSSM) Higgs doublets Hu and Hd, and 
the 16-plet contains the 15 chiral fermions per family of the standard model (SM) 
as well as the right handed neutrino. The implications of this unification have been 
extensively explored over the years [1, 2]. More recently, it has been argued in [3, 4] 
that 50(10) Yukawa unification predicts relatively light (< TeV) gluinos, which can 
be readily tested [5] at the Large Hadron Collider (LHC). The squarks and sleptons 
turn out to have masses in the multi-TeV range. Moreover, it is argued in [3, 4] 
that the lightest neutralino is not a viable cold dark matter candidate, at least in the 
simplest models of 50(10) Yukawa unification. 

Spurred by these developments we have investigated t-b-r Yukawa unification 
[4, 6, 7] in the framework of supersymmetric 5f/(4)e x SU{2)l x SU{2)r [8] (4-2-2, 
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for short). The 4-2-2 structure allows us to consider non- universal gaugino masses 
while retaining Yukawa unification. An important conclusion reached in [4, 6] is 
that with same sign non-universal gaugino soft terms, Yukawa unification in 4-2-2 
is compatible with neutralino dark matter, with gluino co-annihilation [4, 6, 9, 10] 
playing an important role. By considering opposite sign gauginos with /i < 0, M2 < 
0, M3 > (where /z is the bilinear Higgs mixing term, and M2 and M3 are the soft 
supersymmetry breaking gaugino mass terms corresponding respectively to SU{2)l 
and SU{3)c) in [7] we have shown that Yukawa coupling unification consistent with 
known experimental constraints is realized in 4-2-2. With fi < and opposite sign 
gauginos, Yukawa coupling unification is achieved for mo > 300 GeV, as opposed to 
mo ^ 8 TeV for the case of same sign gauginos, by taming the finite corrections to 
the b-quark mass. By considering gauginos with M2 < and M3 > and < 0, we 
can obtain the correct sign for the desired contribution to {g — 2)^ [11]. This enables 
us to simultaneously satisfy the requirements of t-b-r Yukawa unification, neutralino 
dark matter and {g — 2)^, as well as a variety of other known bounds. 

Encouraged by the abundance of solutions and coannihilation channels available 
in the case of Yukawa unified 4-2-2, it is natural to try to further constrain this 
model. One possible way is to impose unification of t-b-r Yukawa couplings with the 
MSSM gauge couplings at Mqut- This is partially inspired from the observation that 
at Mgut, the unified gauge coupling for the MSSM with TeV scale supersymmetry 
is ~ 0.7, while the corresponding third generation Yukawa couplings are of order 
0.6. This suggests that the origin of Yukawa couplings and gauge interaction may 
be closely related, and indeed higher dimensional supersymmetric models have been 
constructed that predict gauge- Yukawa unification (GYU) [12, 13, 14]. We will briefiy 
summarize one such model later in the paper. The phenomenology of this idea was 
studied in [14], where it was shown how a suitable choice of low scale SUSY threshold 
corrections can yield GYU condition in principle, without precisely specifying the 
origin and values for the soft SUSY breaking parameters. 

The main purpose of this paper is to extend the 4-2-2 discussion to the case of 
GYU in the framework of gravity mediated SUSY breaking scenario. In Section 2 we 
briefiy describe the Yukawa unified 4-2-2 model and the boundary conditions for the 
soft supersymmetry breaking (SSB) parameters employed in our scan. In Section 3 
we summarize the scanning procedure and the various experimental constraints that 
we impose. In Section 4 we discuss threshold corrections to the Yukawa couplings 
and summarize from previous studies the findings pertaining to Yukawa unification. 
We also present new results in this section for Yukawa unification for the case of 
> 0, M2 > 0, M3 < in this section. In Section 5 we discuss GYU with MSSM as 
the low energy theory. We first describe a concrete model that breaks to SUSY 4-2-2 
at Mqut and yields gauge-t-6-r Yukawa unification condition. We then proceed to 
discuss the role played by threshold corrections to 6yt in order to obtain GYU. The 
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important role of the top quark mass in implementing GYU is also emphasized. In 
Section 6 we present our results and highlight some of the predictions of the GYU 
4-2-2 model. The correlation between direct and indirect detection of dark matter 
and the gauge- Yukawa unification condition is presented in Section 7 where we also 
display some benchmark points. Our conclusions are summarized in Section 8. 

2 The 4-2-2 model 

In 4-2-2 the 16-plet of 50(10) matter fields consists of i) (4, 2, 1) and V'c (4, 1, 2). The 
third family Yukawa couphng ipci^H, where H{1, 2, 2) denotes the bi-doublet (1,2,2), 
yields the following relation valid at Mqut, 

Yt = Y, = Yr = Y,^. (1) 

In a realistic scenario we can expect corrections to Eq.(l) arising, say, from higher 
dimensional operators. We will assume that these are sufficiently small so that Eq.(l) 
is valid within a few percent or so. 

Supplementing 4-2-2 with a discrete left-right (LR) symmetry [8, 15] (more pre- 
cisely C-parity) [16] reduces the number of independent gauge couplings in 4-2-2 
from three to two. This is because C-parity imposes the gauge coupling unifi- 
cation condition {g^ = Qr) at Mqut- We will assume that due to C-parity the 
SSB mass terms, induced at Mqut through gravity mediated supersymmetry break- 
ing [17] are equal in magnitude for the squarks and sleptons of the three families. 
The tree level asymptotic MSSM gaugino SSB masses, on the other hand, can be 
non-universal from the following consideration. From C-parity, we can expect that 
the gaugino masses at Mqut associated with SU{2)l and SU{2)r are the same 
(M2 = M2 = M2). However, the asymptotic S'f/(4)c and consequently SU {3)c 
gaugino SSB masses can be different. With the hypercharge generator in 4-2-2 given 
by Y = a/2/5 {B — L) + a/3/5 I3R, where B — L and I^r are the diagonal generators 
of SU{4:)c and SU{2)r, we have the following asymptotic relation between the three 
MSSM gaugino SSB masses: 

Ml = + ^M3. (2) 

The supersymmetric 4-2-2 model with C-parity thus has two independent param- 
eters (M2 and M3) in the gaugino sector. In order to implement Yukawa unification it 
turns out that the SSB Higgs mass terms must be non-universal at Mqut- Namely, 
< fnjj^ at Mqut, where mH^^mHj is the up (down) type SSB Higgs mass 
term. Phenomenological studies of the MSSM typically resort to 'just-so' splitting 
(see Blazek, Dermisek and Raby in [2]) for the MSSM Higgs doublets, while remark- 
ing in passing that such a splitting may arise, for example, from D terms. D terms. 
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however, induce splitting in the squarks and sleptons as well. It is possible to imagine 
a simple mechanism to implement just-so Higgs SSB mass splitting in either the 4-2-2 
or SO{10) models. 

We need Higgs fields other than H{1, 2, 2) in order to complete the model. For in- 
stance, $(4, 1, 2) + j>(4, 2, 1) field may be used to break 4-2-2 to the SM gauge group. 
In the 4-2-2 model, just-so Higgs splitting may be understood by writing an SU{2)ji 
violating bilinear term between the up and down type Higgs doublets. This may be 
done, for instance, by considering a super-heavy Higgs A that transforms as (1, 1, 3) 
under 4-2-2. If the MSSM doublet comes from H{1, 2, 2) under 4-2-2, the superpoten- 
tial will have the interaction HHA. We may consider the SSB trilinear term A/^HHA 
and assume that A has a non-zero VEV, (A) = Vdiag{l, —1) with V ~ (O) TeV. 
We can thus achieve the desired splitting in the MSSM Higgs doublets. Likewise, we 
can explain just-so Higgs splitting in 5*0(10) by using the 45-dimensional Higgs since 
it contains the A(l, 1, 3) of 4-2-2. We may consider the non-renormalizable coupling 
10 ■ 10 ■ 45^ (10 •10-45 is not allowed as 45 is a two-index anti-symmetric represen- 
tation) in the superpotential, where the 10 is the 10-dimensional representation that 
has the MSSM Higgs doublets. We may get just-so Higgs splitting by writing the 
corresponding SSB term with a suitable choice of Yukawa coupling or VEV for the 
45. 

The fundamental parameters of the 4-2-2 model that we consider are as follows: 

mo, mH^,mHa, M2, M3, Aq, tan/3, sign(/i). (3) 

Here mo is the universal SSB mass for MSSM sfermions, Aq is the universal SSB 
trilinear scalar interaction (with the corresponding Yukawa coupling factored out), 
tan 13 is the ratio of the vacuum expectation values (VEVs) of the two MSSM Higgs 
doublets, and the magnitude of fi, but not its sign, is determined by the radiative 
electroweak breaking (REWSB) condition. Although not required, we will assume 
that the gauge coupling unification condition gs = gi = g2 holds at Mqut in 4-2-2. 
Such a scenario can arise, for example, from a higher dimensional S'O(IO) [18] or 
SU (8) [12] model after suitable compactification. 

3 Phenomenological constraints and scanning pro- 
cedure 

We employ the ISAJET 7.80 package [19] to perform random scans over the parameter 
space listed in Eq.(3). In this package, the weak scale values of gauge and third 
generation Yukawa couplings are evolved to Mqut via the MSSM renormalization 
group equations (RGEs) in the DR regularization scheme. We do not strictly enforce 
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the unification condition = gi = g2 at Mgut? since a few percent deviation from 
unification can be assigned to unknown GUT-scale tliresliold corrections [20]. Tfie 
difference between gi{= (72) and g^ at Mqut is no worse tlian 4%. If neutrinos acquire 
mass via Type I seesaw, tlie impact of tlie neutrino Dirac Yukawa coupling on tlie 
RGEs of tfie SSB terms, gauge couplings and the third generation Yukawa couplings 
is significant only for relatively large values (~ 2 or so). In the GYU 4-2-2 model 
we expect the largest (third family) Dirac Yukawa coupling to be comparable to the 
gauge couplings (~ 0.6 at Mqut)- Therefore, we do not include the Dirac neutrino 
Yukawa coupling in the RGEs. 

The various boundary conditions are imposed at Mqut and all the SSB parame- 
ters, along with the gauge and Yukawa couplings, are evolved back to the weak scale 
Ml- In the evaluation of Yukawa couplings the SUSY threshold corrections [21] are 
taken into account at the common scale Msusy = ^m^^m^^. The entire parameter 
set is iteratively run between Mz and Mgut using the full 2-loop RGEs until a stable 
solution is obtained. To better account for leading-log corrections, one-loop step-beta 
functions are adopted for gauge and Yukawa couplings, and the SSB parameters rrii 
are extracted from RGEs at multiple scales m, = mi{mi). The RGE-improved 1- 
loop effective potential is minimized at an optimized scale Msusy, which effectively 
accounts for the leading 2-loop corrections. Full 1-loop radiative corrections are in- 
corporated for all sparticle masses. 

The requirement of REWSB [22] puts an important theoretical constraint on the 
parameter space. Another important constraint comes from limits on the cosmologi- 
cal abundance of stable charged particles [23] . This excludes regions in the parameter 
space where charged SUSY particles, such as fi or ti, become the lightest supersym- 
metric particle (LSP). We accept only those solutions for which one of the neutralinos 
is the LSP and saturates the WMAP (Wilkinson Microwave Anisotropy Probe) dark 
matter rehc abundance bound. 

We have performed random scans for the following parameter range: 

< mo, niH^ , niH^ < 20 TeV 
-2TeV < M2 < 2 TeV 
-2TeV < M3 < 2 TeV 
45 < tan /3 < 55 
-3 < Aq/tjiq < 3 

/i<0,;U>0 (4) 

where = 173.3 ± 1.1 GeV [24] is the top quark pole mass. The value of the top 
quark mass is very crucial, as we shall see later, for GYU. We use mb(m^) = 2.83 
GeV which is hard-coded into ISAJET. The above choice of parameters is influenced 
by our previous experience with the 4-2-2 model. 
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In scanning the parameter space, we employ the Metropohs-Hastings algorithm as 
described in [25]. All of the collected data points satisfy the requirement of REWSB, 
with the neutralino in each case being the LSP. We direct the Metropolis-Hastings 
algorithm to search for solutions with GYU. After collecting the data, we impose the 
mass bounds on all the particles [23] and use the IsaTools package [26] to implement 
the following phenomenological constraints on points that have GYU to within 20%: 



nih (lightest Higgs mass) 
BR{Bs /i+/i-) 
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[30] 
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In the case of A{g — 2)^, we only require that the GYU 4-2-2 model does no worse 
than the SM. However, we do give examples of solutions that satisfy the A{g — 2)^/2 
constraint to within 3a. 



4 Threshold corrections and Yukawa unification 

The SUSY threshold corrections to the top, bottom and tau Yukawa couplings play 
a crucial role in t-b-r Yukawa coupling unification. In general, the bottom Yukawa 
coupling yb can receive large threshold corrections, while the threshold corrections to 
Ut are typically smaller [21]. The scale at which Yukawa coupling unification occurs 
is set equal to Mqut, the scale of gauge coupling unification. Consider first the 
case ?/t(MGUT) ~ 2/r(^GUT)- The SUSY correction to the tau lepton mass is given by 
SitLt- = V cos (36yr- For the large tan /3 values of interest here, there is sufficient freedom 
in the choice of Syr to achieve yt ~ Vr at Mqut- This freedom stems from the fact 
that cos/5 ~ 1/ tan/3 for large tan/3, and so we may choose an appropriate Syr and 
tan/3 to give us both the correct r lepton mass and yt ~ yr- The SUSY contribution 
to Syb has to be carefully monitored in order to achieve Yukawa coupling unification 

We choose the sign of Syi {i = t, b, r) from the perspective of evolving yi from 
Mgut to Mz- With this choice, Syb must receive a negative contribution (—0.27 < 
Syb/Ub ^ —0.15) in order to reahze Yukawa coupling unification [7]. This is a narrow 
interval considering the full range of —0.4 < Syb/yb ^ 0.6 that we found in the data 
that we collected. The dominant contribution to Syb comes from the finite corrections 
of the gluino and chargino loops, and in our sign convention, is approximately given 
by [21] 
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finite 



(5) 



127r2 m? 327r2 m? 
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where is the strong gauge couphng, rrig is the gluino mass, mi and are the 
hghter sbottom and stop masses, and At is the top trihnear (scalar) couphng. 

The logarithmic corrections to i/b are positive, which leaves the finite corrections 
to provide for the correct overall negative Syb in order to realize Yukawa unification. 
The gluino contribution (Eq.(5)) is positive for > and same sign gaugino soft 
mass terms. Thus, the chargino contribution (Eq.(5)) must play an essential role to 
provide the required negative contribution to 6yb. This can be achieved with suitably 
large values of both mo and At. This large value of mo and At is the reason behind 
the requirement of mo > 6TeV and Ao/mo ~ —2.6 in the 5*0(10) model discussed 
in [3]. The parameter tan/3 also lies in a narrow range 48 < tan /3 < 52. A similar 
trend was shown in [4] for the 4-2-2 model with same sign gauginos and n > 0. The 
latter model displays Yukawa coupling unification consistent with WMAP data via 
bino-gluino coannihilation. 

In an S'O(10)-like [3] model with same sign gauginos, the case /i < is not 
favored because of the negative contribution to A{g — 2)fj, oc fiM2 which, instead, 
needs to be positive. Therefore, while the Yukawa unified SO{10) and 4-2-2 [4] (with 
/i, M2, M3 > 0) models do not provide the required contribution to A{g — 2)^ because 
of heavy sparticles, they do no worse than the SM in this respect. 

One can improve the situation immensely by considering the case of opposite 
sign gaugino soft terms which is allowed by the 4-2-2 model. We showed in [7] the 
parameter space corresponding to /x, M2 < 0, M3 > that gives Yukawa coupling 
unification with a sub-TeV sparticle spectrum which is consistent with all known 
experimental bounds including A{g — 2)^. Another possibility is to consider /i, M2 > 
0, M3 < (the parameter space for this case has not been previously discussed in the 
literature). This becomes possible because the gluino contribution to 6yb is of the 
correct (negative) sign. 

In order to quantify Yukawa coupling unification, following [3], we define the 
quantity R as. 



Thus, i? is a useful indicator for Yukawa unification with R < 1.1, for instance, 
corresponding to Yukawa unification to within 10%, while R = 1.0 denotes 'perfect' 
Yukawa unification. 



In Figure 1 we show a comparison between the four model types considered, i. e. 
the 50(10) model and 4-2-2 models with {/i > 0, M2 > 0,M3 > 0}, {fi < 0,M2 < 



R 



max{yt,yb,yT) 
^MyuVb^yr) 



(6) 
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0, Ms > 0} and {fi > 0, M2 > 0, M3 < 0}. We show plots in the mo - ^oMo plane 
for these models. Gray points are consistent with REWSB and Xi LSP. Orange, light 
blue and purple points are subsets of gray points with R < 1.2, 1.1, 1.02 respectively. 
As previously explained, in the SO (10) and 4-2-2 models with > and same sign 
gauginos, Yukawa coupling unification can only be achieved for large values of rriQ. 
Also, the value of Ag/mo is very restricted. On the other hand, with opposite sign 
gauginos one can realize very credible Yukawa unification with relatively small tuq 
values. Also as previously described, since the gluino loop provides the required 6yb, 
Aq/itlq is no longer restricted by Yukawa coupling unification and can vary over a 
very wide range. 

In Figure 2 we display an interesting difference between the same sign and opposite 



(a) SO{lQ) (b) Ai > 0, M2 > 0, M3 > 




-3-2-10123 -3-2-10123 

Ao/ mo Aol mo 



Figure 1: Plots in the mo - Aq/tjiq plane for the 5*0(10) model and three classes of 
4-2-2 models. Gray points are consistent with REWSB and x? LSP. Orange, light 
blue and purple points are subsets of gray points with R < 1.2, 1.1, 1.02 respectively. 
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sign gaugino cases in the IM3I - mo plane. Shown in gray are points that satisfy the 
requirements of REWSB and Xi LSP. In green, blue and orange, we show points that 
further satisfy Yukawa unification to within 20%, 15% and 10% respectively. Red 
points satisfy the particle mass bounds in addition to having R < 1.1. The trend of 
a lower M3 in the case of same sign gauginos is very apparent if we require Yukawa 
coupling unification. The reason for this is again that we need to suppress the finite 
correction to the bottom quark mass coming from the gluino loop (see Eq.(5)). The 
case of opposite sign gauginos, in stark contrast, shows that essentially any value of 
IM3I is acceptable as far as Yukawa coupling unification is concerned. The orange 
region in the bottom left in this case is excluded because of the lower bound on the 
gluino mass. It is also instructive to consider the rrig - At plane in the case of same 
sign gauginos shown in Figure 3. In this figure, rrig is the physical gluino mass and 
At is the value of the top trilinear (scalar) coupling at the scale Q = ^mi^mij^. The 
color coding is the same as in Figure 2. This figure shows that a lighter gluino is 
required for Yukawa coupling unification with a smaller absolute value of At. This 
again reflects the fact that for Yukawa coupling unification, we need to suppress the 
gluino contribution to dyi, in favor of the chargino contribution. 

Finally, it is interesting to note that with opposite sign gauginos, the MSSM 
parameter tan (5 varies over a wider range, 44 < tan /3 < 54. With same sign gauginos, 
on the other hand, 48 < tan (5 < 52. Among other interesting features of a Yukawa 
unified model with opposite sign gauginos is the presence of various channels for 



(a) Same Sign Gauginos (b) Opposite Sign Gauginos 




' '■ ■ ■ — ■■ T . II .1 I — OS Q r. . . . L f.":T' I ... -1- ... I - ...I 
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Figure 2: Plots in the IM3I - tuq plane for same sign (left) and opposite sign (right) 
gaugino (including points for both {/i < 0, M2 < 0, M3 > 0} and {/i > 0, M2 > 
0, M3 < 0}) cases. Gray points are consistent with REWSB and Xi LSP. Green, blue 
and orange points are subsets of gray points with R < 1.2, 1.15, 1.1 respectively. Red 
points satisfy particle mass bounds in addition to i? < 1.1. 
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A(TeV) 



Figure 3: Plot in the rrig - At plane plane for same sign gauginos. Color coding same 
as in Figure 2. 



realizing the desired Xi ^^^lic density. In particular, in [7] we showed the existence 
of stau coannihilation, bino-wino coannihilation, gluino coannihilation and CP-odd 
Higgs resonance solutions for the case (/i < 0, M2 < 0, M3 > 0). In contrast, for the 
Yukawa unified 5*0(10) model, only the light Higgs resonance solution is consistent 
with the WMAP relic density. It is interesting to note that in the case of 4-2-2 
models (with /i > 0, M2 > 0, M3 > 0), it is not possible to get the well-known 
stau (f) coannihilation channel. This is because in the f mass^ matrix the diagonal 
terms are proportional to mg, whereas the off-diagonal terms are proportional to 
At-itLt-, where Aj. is the low-scale value of the tau trilinear (scalar) coupling. For 
/i > 0, M2 > 0, M3 > 0, one needs a heavy mo in order to realize Yukawa couphng 
unification. One therefore needs A^- ~ m^/rrir in order for the off-diagonal terms to 
contribute to give a small stau mass mf ~ ""^^o, where Xi is the lightest neutralino. 
This, clearly, is not possible for large uiq values. The parameter space of the 5*0(10) 
model is a subset of the 4-2-2 model with {/i > 0, M2 > 0, M3 > 0} and so these 
remarks apply to 50(10) as well. 



5 Gauge- Yukawa unification 

In this section, we first describe a specific model where GYU may happen. We then 
move on to discuss SUSY particle thresholds and their effects on analyzing GYU. It 
is helpful to define, in analogy with R, a parameter GY that quantifies GYU; 

= max(gfi,gf2,g'3,yf,yb,Z/r) 
min(^i,^2,^3,2/t,2/6,yr) 



10 



5.1 Model for Gauge- Yukawa unification 

A six dimensional model realizing unification of the gauge couplings {gi, g2, gs) 
and the third family Yukawa couplings {yt, yb, yr) was presented in [12]. It has 
SU{8) gauge symmetry with N=2 SUSY, which corresponds to N=4 SUSY in 4D, 
and thus only the gauge multiplet can be introduced in the bulk. The 6D N=2 
gauge multiplet, expressed in terms of 4D, N=4 gauge multiplet, contains the vec- 
tor multiplet V{An, A) and three chiral multiplets in the adjoint (63-dimensional) 
representation of the gauge group. The 63-dimensional gauge multiplet contains 
the gauge bosons (and their superpartners), while the three 63-dimensional chi- 
ral multiplets contain the third family matter fermions and the Higgs bosons plus 
their superpartners. The two extra dimensions are compactified on the orbifold 
T^/Zg, and a suitable choice of the Zq transformation matrix breaks SU{8) down 
to SU{A) X SU{2)l X SU{2)r x U{lf. The theory reduces to 4D N=l SUSY 4-2-2 
model and two additional U{1) factors. The massless modes after compactification 
are the 4-2-2 gauge fields, (15, 1, 1), (1, 3, 1), (1, 1, 3) two singlet vector fields (1, 1, 1) 
and (1, 1, 1), third-family matter fermions = (4, 2, 1)2 q and = (4, 1, 2)_2 
and the bi-doublet Higgs fields. Hi = (1, 2, 2)^ ^ and H2 = (1, 2, 2)^ 
The trilinear coupling for the chiral multiplets 

S = j(fx J (f9 2Ti(^-V2gG^<^,<!?']^ +h.c. (8) 

includes the third family Yukawa interaction terms 

S = j d^x j (fdyfs^LHi-^R + h.c. (9) 

In Eq. (8), E, $, ^'^ are chiral multiplets containing the third family chiral fields, ^ l 
and and the bi-doublet Higgs fields. Hi and H2, and g^ and yQ are the 6D gauge 
and Yukawa couplings. Eqs. (8) and (9) lead to g^ = yQ with proper renormalization 
of the kinetic terms. Integrating out the two extra dimensions, we obtain y^ = g4 for 
the 4D coupling leading to 

9i = 92 = g3 = yt = yb = yA= y^^n (10) 

at the compactification scale (Mc) which we identify with the four dimensional uni- 
fication scale Mgut- We assume that the 4-2-2 symmetry, as well as the two extra 
U (1) are broken at Mqut to the SU{3)c x SU {2)^x11 (l)y using suitable Higgs vevs 
on the brane. We further assume that the breaking of the two extra U{1) symmetries 
does not cause D term splittings so that the universality of squark mass^ terms at 
Mgut is preserved. The first and second families are treated as brane fields to cancel 
the brane localized gauge anomalies. The Yukawa couplings for the first and second 
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families are suppressed by a large volume factor, but there is no good reason as to 
why the mass of the first family is hierarchically small. The particle spectrum below 
Mgut is the same as in MSSM. For related discussions see Ref. [13]. 

5.2 SUSY thresholds and Gauge- Yukawa unification 

As previously mentioned, the bottom quark and tau lepton Yukawa couplings receive 
larger threshold corrections than yt. Since the gauge coupling is more or less fixed 
~ 0.69, and since Yukawa coupling unification typically occurs for ^ ^ ~ 0.6, 
the quantity GY ~ 1.15 for the Yukawa unified models discussed above. If we desire 
to impose GYU on our models, yt becomes the bottleneck for a given top mass as 
typically and y^. can be made larger than yt by a suitable choice of the SUSY 
spectrum. In particular, the values of yb and yr at Mgut can be pushed up to 
?/b(MGUT) ~ 1-25'GUT and i/r(MGUT) ~ 1-2(?gut, where (^gut is the value of the 
unified gauge coupling at Mgut- The leading SUSY threshold correction to the top 
quark mass is given by [21] 

x„ finite^ 9l f^m-g tan 13 
^ 127r2 m? 

In our sign convention (evolving the couplings from Mgut to M^), a negative con- 
tribution to 5yt is preferred. Naively, a larger negative contribution allows for a larger 
yt{MG\ji). However, in the case of same sign gauginos with /i > 0, we get a positive 
contribution to dyt-, in which case a large mo value is required. The requirement of a 
large can be argued from just requiring Yukawa coupling unification. The signifi- 
cance of looking at the sign of the correction to 5yt in this case is the realization that 

1.3 
1.2 

1.1 
1 

S 0.9 
=S 0.8 

0.7 

0.6 

0.5 

0.4 



Figure 4: Gauge- Yukawa unification to within 6% in SU{4)c x SU{2)l x SU{2)r 
(4-2-2). 
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it may not be possible to achieve (more or less) gauge- Yukawa unification at all. We 
see that GY > 1.13 in the data that we have collected. In the case of opposite sign 
gauginos, on the other hand, our choice of the sign of /x gives a negative contribution 
to 6yt. We should, therefore, expect that GYU is allowed in the case of opposite sign 
gauginos. 

In Figure 4 we show the evolution of the gauge couplings and the third generation 
Yukawa couplings that unify to within 6% in the 4-2-2 model. The spectrum for this 
point is given as Point 4 in Table 1. 

5.3 Gauge- Yukawa unification and rrit 

It is perhaps not too surprising that the parameter GY, a measure of GYU, depends 
sensitively on the top quark mass rrit. It is, therefore, instructive to study how 
GYU is affected as one varies rrit. We plot in Figure 5 GY as a function of rrit. As 
expected, GYU prefers a larger top mass, with near perfect unification possible for 
rrit = 177GeV. We next discuss GYU allowing for a la variation in the top mass. 

6 Gauge- Yukawa unification and sparticle spectro- 
scopy 

We present here the results of the scan over the parameter space listed in Eq.(4) after 
allowing for rut to vary within la of its central value. In Figure 6 we show results 
in the M3 - M2 plane. As previously explained, gauge- Yukawa unification prefers 
opposite sign gauginos. We emphasize this by only showing the two cases of relative 
sign gauginos and leaving two quadrants empty in Figure 6. (Same sign gauginos GY 



1.2 
1.18 




171 172 173 174 175 176 177 

?77,(GeV) 

Figure 5: Plot of GY versus rut. The vertical lines correspond to nit = 172.2, 173.3, 
174.4 and 175.5 GeV. 
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unification of order 10% or higher. See later.) The gray points are consistent with 
REWSB and Xi LSP, while the green points also satisfy the particle mass bounds and 
constraints from BR{Bs — ?■ ^'^^^), BR{h — )■ 57) and BR{Bu — ?■ tVt). In addition, we 
require that the green points fare no worse than the SM as far as — is concerned. 
The blue points belong to the subset of green points that satisfies the WMAP bounds 
on Xi dark matter abundance. Points in red represent the subset of blue points that 
satisfies gauge- Yukawa coupling unification to within 10%. We also show the lines 
M3 = — 6.3M2 and M3 = — O.I2M2. The slopes of these lines indicate bino-wino 
and bino-gluino coannihilation in the M3 - M2 plane. If we start off with a universal 
gaugino mass at Mgut, we get M2/Mi(Q) ^ ±1.89 and M^/Mi^Q) ^ ±4.67, where 
the negative sign is for the case of opposite sign gauginos. Therefore, in order to get 
bino-wino coannihilation we should set Mi/M2(Mgut) ~ ±1.89. Substituting this 
ratio of Mi and M2 in Eq. (2) we can infer that bino-wino coannihilation will occur 
for M3 ^ — 6.3M2. A similar calculation shows that for bino-gluino coannihilation 
we should set M3 ^ -O.I2M2. 

In Figure 7 we show the results in the GY - ttiq and GY - tan /3 planes. There 
is no visible distinction between the two classes of opposite sign gaugino models 
in these two planes, which is why we plot the data from the two sets in the same 



Figure 6: Results in the M3 - M2 plane. Gray points are consistent with REWSB and 
Xi LSP. Green points satisfy particle mass bounds and constraints from BR{Bs — )■ 
fi^fi^), BR{h —7- 57) and BR{Bu — )■ tz/^). In addition, we require that green points 
do no worse than the SM in terms of {g — 2)^. Blue points belong to a subset of green 
points and satisfy the WMAP bounds on x^ dark matter abundance. Points in red 
represent the subset of blue points satisfying gauge- Yukawa coupling unification to 
within 10%. We also show the lines 2M3 = -I3M2 and 4IM3 = -6M2 discussed in 
the text. 
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Figure 7: Plots in the GY - mo and GY - tan /3 planes. The two classes of opposite 
sign gaugino models are shown together. Color coding same as in Figure 6. 



figure. The color coding is the same as in Figure 6. It can be seen that a relatively 
large mo(~ 500 GeV) is required even without imposing any of the experimental 
constraints. After including the experimental constraints, we are forced to have 
rriQ ~ 1.5 TeV. This is to be contrasted with the situation depicted in Figure 1 
where tuq ~ 300 GeV suffices for Yukawa coupling unification compatible with all 
constraints. This may be understood from the fact that keeping all other parameters 
fixed, a larger ttiq value tends to push up the value of ?/t(MGUT) closer to the unified 
gauge coupling. Likewise, we must have 46 < tan /3 ^ 54 with an even narrower range 
(47 < tan /3 < 52) if we consider the experimental constraints. Yukawa coupling 



(a) /Lt < 0, M2 < 0, M3 > 
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(b) /X > 0, M2 > 0, M3 < 
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Figure 8: Plots in the GY - Ao/mo plane for the two classes of opposite sign gaugino 
models. Color coding same as in Figure 6. 
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unification, on the other hand, allows for 44 < tan /3 ^ 54 for Yukawa unification 
consistent with experimental constraints. 

In Figure 8 we show plots in the GY - Ao/nio plane for the two classes of opposite 
sign gaugino models. The color coding is the same as in Figure 6. It is evident 
that GYU with /i > prefers Ao/mo < 0, and vice versa. This is different from 
just Yukawa unified 4-2-2 as seen clearly from Figure 1. This stems from the finite 
chargino contribution to Syt, which is proportional to fiAt. In the case of Yukawa 




0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 

771^0 (TeV) 777^0 (TeV) 



Figure 9: Plots in the m^± - m^o, rrig - m^o, ruf - m^o and tua - m^o planes. The gray 
points satisfy the requirements of REWSB and Xi LSP. The blue points, in addition, 
satisfy particle mass bounds and constraints from BR[Bs — )■ fi~^fi~), BR{Bu — )■ rur) 
and BR{b — )■ 57). In addition, we require that these points do no worse than the 
SM in terms of the {g — 2)^ prediction. The red points correspond to GYU to within 
10% in addition to these constraints. We show in the m-± - m^o, irtf, - m^o and - 
m^o planes the unit slope lines that indicate the respective coannihilation channels. 
In the niA - ttl^o plane we show the line mA = 2m^o that signifies the A resonance 
channel. 
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unification, one can have a small mo value for which the chargino contribution is 
sub- dominant. In GYU on the other hand, mg is large as previously explained. This, 
coupled with the fact that we need the threshold correction to 6yb to be negative, 
shows that ^Aq/itiq < is preferred for GYU. 

In Figure 9 we show the relic density channels consistent with GYU in the m-± 
- ^j^o, rrig - m^o, mf - m^o and - m^o planes. The gray points in this figure 
satisfy the requirements of REWSB and x? LSP. The blue points, in addition, satisfy 
the particle mass bounds and constraints from BR{Bs — )■ fi~^fi~), BR{Bu — > rz/,-) 
and BR{b — )■ 57). In addition, we require that these points do no worse than the 
SM in terms of the {g — 2)^ prediction. The red points correspond to GYU to 
within 10% in addition to these constraints. We can see in Figure 9 that a variety 
of coannihilation and annihilation scenarios are compatible with Yukawa unification 
and neutralino dark matter. Included in the - m^o plane is the line = 2m^o 
which indicates that the A funnel region is compatible with GYU. In the remaining 
planes in Figure 9, we draw the unit slope line which indicates the presence of gluino 
and stau coannihilation and bino-higgsino mixed dark matter scenarios. 

Recent results from ATLAS [32, 33] and CMS [31] naively put very stringent 
limits on the gluinop mass of nig > 500 GeV. However, it is shown explicitly in [34] 
that this limit does not apply in general, and specifically, does not apply in the case 
of heavy > TeV squarks. In the gluino coannihlation channel shown in Figure 9 the 
squarks are heavy and are not yet excluded by ATLAS/CMS. The recent results do 
seem to suggest that the gluino coannihilation scenario will soon be tested. 

7 Gauge- Yukawa unification and dark matter de- 
tection 

In light of the recent results by the CDMS-II [35] and XenonlOO [36] experiments, it 
is important to see if GYU, within the framework presented in this paper, is testable 
from the perspective of direct and indirect detection experiments. The question of 
interest is whether /i ~ Mi is consistent with GYU, as this is the requirement to 
get a bino-higgsino admixture for the lightest neutralino which, in turn, enhances 
both the spin dependent and spin independent neutralino-nucleon scattering cross 
sections [37]. In Figure 10 we show the spin independent and spin dependent cross 
sections as a function of the neutralino mass. In the case of spin independent cross 
section, we also show the current experimental bounds and future reach of the CDMS 
and Xenon experiments. The color coding is the same as in Figure 9. A small region of 
the parameter space consistent with GYU and the experimental constraints discussed 
in Section 3 (red points in the figure) is at the exclusion limits set by the current 
CDMS and XENON experiments. Thus, the ongoing and planned direct detection 
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?7i^o(GeV) 771^0 (GeV) 

Figure 10: Plots in the cxsi - m^o and ctsd- pl^-nes. Color coding is the same as in 
Figure 9. In the asi - m^o plane we show the current bounds (solid lines) and future 
reaches (dotted lines) of the CDMS (black lines) and Xenon (red lines) experiments. 
In the (TsD - f^x^ plane we show the current bounds from Super K (black line) and 
IceCube (dotted black line) and future reach of IceCuce DeepCore (dotted red hne). 



experiments will play a vital role in testing GYU models. 

In the case of spin dependent cross section, we show in Figure 10 the current 
bounds from the Super-K [38] and IceCube [39] experiments and the projected reach 
of IceCube DeepCore. The current Super-K and IceCube bounds are not stringent 
enough to rule out anything. However, from Figure 10 we see that the future IceCube 
DeepCore experiment will be able to constrain a significant region of the parameter 
space. 

In Table 1 we present some benchmark points for the 4-2-2 GYU model. All of 
these points are consistent with neutralino dark matter and the constraints mentioned 
in Section 3. Point 1 represents the best GYU that we have found and corresponds to 
the A funnel region. Points 2 and 3 correspond to the gluino and stau coannihilation 
channels, while for Point 4 bino-Higgsino mixing plays a major role in giving the 
correct dark matter relic density. As expected, both the spin independent and spin 
dependent cross sections of the neutralinos on protons are larger for Point 4. Note 
that Point 3 also satisfies the lower bound on A (^f — 2)^. Finally, point 5 represents 
the best GYU solution that we found in the case of same sign gauginos [GY ~ 1.14). 
All of the points shown in this Table are currently allowed by ATLAS/CMS [34]. 

8 Conclusions 

Guage- Yukawa unification (GYU) at Mqut, a natural extension of four dimensional 
gauge unification, is implemented using higher dimensional theories in which the 
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Point 1 


Point 2 


Point 3 


Point 4 


Point 5 


mo 


2063 


3246 
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1769 


7171 


Ml 


747 


1034 


-418 


985 


583 


M2 


1742 


1819 


-1455 


1938 


939 


M3 


-744 


-143 


1138 


-443 


49 


tan /3 


50 


48 


47 


50 


53 


Ao/mo 


-1.73 


1.94 


-0.19 


-0.61 


-2.53 


rriHu 


2191 


1162 


657 


1328 


4557 


rriHd 


2797 


3286 


1294 


2330 


6722 


rrit 


174.3 


174.1 


174.2 


174.4 


173.1 


sgn /X 


+1 


+1 


-1 


+1 


+1 


rrih 


117 


119 


119 


116 


121 


niH 




1 Ton 

1 /oy 


694 


1 1 no 
iiUz 


HOT 


rriA 


594 


1728 


689 


1095 


r\ d 

983 


mH± 


DUO 




( Ui 


i iUO 


OOQ 


m^o 

Al,2 

A3, 4 

m-± 

Al,2 


340, 528 


459, 1530 


193, 1040 


428, 492 


296, 907 


^9Q 14Q9 




in^R i9^(S 

lUtJO, XZitJU 


^{V\ l(^9Q 

OUO, lUZc/ 




534, 1478 


1531, 2709 


1050, 1247 


497, 1618 


r\r\r\ r^r^c^c 

909, 6686 




1750 


516 


2503 


1128 


340 




2732, 2485 


3445, 3209 


2451, 2260 


2319, 1946 


7186, 7110 




1355, 1793 


1788, 2091 


1846, 2091 


1153, 1643 


1948, 2607 


mr 

0l,2 


2733, 2511 
1336, 1781 


3445, 3276 
1485, 2065 


2452, 2272 
180L 2074 


2321, 1982 
924, 1635 


7186, 7195 
2407, 2852 




2335 


3412 


1177 


2148 


7161 




1841 


2814 


1048 


1835 


5218 




2336, 2115 


3412, 3335 


1181, 784 


2149, 1854 


7160, 7251 




540, 1836 


1911, 2815 


202, 1059 


947, 1833 


2129, 5204 


(757 (pb) 


9.1 X 10-*^ 


4.7 X 10-^2 


2.7 X 10-1° 


2.5 X 10-^ 


1.1 X 10-12 


(pb) 


5.6 X 10-6 


5.4 X 10-^° 


9.0 X 10-^ 


3.4 X 10-^ 


7.8 X 10-12 


0.09 


0.1 


0.11 


0.08 


0.10 


R 


1.05 


1.07 


1.08 


1.04 


1.13 


GY 


1.05 


1.09 


1.09 


1.06 


1.14 



Table 1: Point 1 is the best GYU we found corresponding to the A funnel region. 
Points 2 and 3 respectively correspond to the gluino and stau coannihilation channels, 
while for Point 4 the LSP is a bino-Higgsino admixture. Point 5 represents the best 
GYU solution we found with same sign gauginos and corresponds to gluino NLSP. 

gauge and third {t-b-r) family matter supermultiplets are unified. One of the simplest 
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realizations of this idea gives rise, after compactification, to the well-known symmetry 
group SU{4:)c X SU{2)l x SU{2)r. GYU in this framework strongly prefers gaug- 
ino masses M2 and M3 with opposite signs, and it also shows some preference for a 
top mass that is slightly higher than 173.3 GeV, its current central value. We have 
explored the fundamental parameter space of GYU models and identify a number 
of benchmark points that are compatible with a large variety of experimental con- 
straints, including the WMAP bound on neutralino dark matter and {g — 2)^. One 
of the more intriguing GYU compatible solutions corresponds to the gluino NLSP 
scenario which can be tested at the LHC. 
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